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_ ■ I. INTRODUCTION 

> ■ 

' Impulsive pp- waves (plane fronted gravitational waves with parallel rays) can be described by a metric of the form 
O . 

O ■ ds^ = 6{u) f{x,y)du'^ - dudv + dx^ + dy^ , (1) 

\0 ; 

, where {u, v) and {x, y) is a pair of null and (transverse) Cartesian coordinates respectively, and / denotes the profile 
function subject to the field equations. Hence the spacetime is flat everywhere except for the null hypersurface 
— ^ , u = 0, where it has a (5-like pulse modelling a gravitational shock wave. Such geometries arise most naturally as 
O ' ultrarelativistic limits of boosted black hole spacetimes of the Kerr-Newman family (as shown by various authors |2|-|j] ) 
and multipole solutions of the Weyl family ||^. Also they play an important role in particle scattering at the Planck 
?H . scale (see [|| and references therein). 

^f^' There have been intrinsic descriptions of impulsive pp-waves, viz. by Penrose and by Dray and t'Hooft 0, which 
^ , essentially consist in glueing together two copies of Minkowski spacetime with a warp across the null hypersurface 
' u — 0. Penrose also introduced a different coordinate system in which the components of the metric tensor are actually 
rS continuous. However, the transformation relating the coordinates used in (^ to these new ones is discontinuous (for 
the general form of the transformation see Q) and therefore -strictly speaking- the differential structure of the 
manifold is changed. In this paper we stick to the original distributional form of the metric, motivated by the fact 
that physically, i.e. in the ultrarelativistic limit, the spacetime arises that way (cf. the approaches of PJTc|]). For recent 
work on pp-waves using the continuous form of the metric see pT| . 

We describe the geometry of impulsive pp-waves entirely in the distributional picture using the framework of 
Colombeau's generalized functions, thereby generalizing previous work p^ . As discussed there in detail, the geodesic 
as well as the geodesic deviation equation for impulsive pp-waves involve formally ill-defined products of distributions, 
due to the nonlinearity of the equations and the presence of the Dirac (5-function in the space time metric. However, as 
was also shown in , one can overcome these difficulties using a careful regularization procedure which, while math- 
ematically sound, corresponds to the physical idea of viewing the impulsive wave as the limiting case of a sandwich 
wave of ever decreasing support but constant (integrated) strength. More precisely, regularizing the (S-distribution 
by a "model S-net" (i.e., a net Pe{x) := p{xe~^), where p is a smooth function with support contained in the 
interval [—1, 1] satisfying J p = 1), it was shown that the solutions to the smoothened equations possess regularization 
independent weak limits. These distributional "solutions" fit perfectly into the physically expected picture showing 
that the geometry of impulsive pp-waves can be described consistently using the distributional form of the metric. 
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The reliability of the results is guaranteed by making use of regularization techniques instead of introducing "multi- 
plication rules" into Schwartz linear disribution theory (cf. the discussion at the end of Sec. 2 in or |Q for general 
remarks) . 

However, the "solutions" obtained by this naive regularization procedure exhibit a mathematically highly unsatisfac- 
tory feature. They do not obey the original distributional equations (unless, again, one is willing to impose certain 
"multiplication rules"), as is common to such situations. Hence -strictly speaking- this approach does not provide 
a reasonable solution concept for the equations under consideration. Such a notion is available in the nonlinear the- 
ory of generalized functions [p^-p^ due to J. F. Colombeau, where one has -loosely speaking- a rigorous system 
of bookkeeping on the regularizing sequences. Recently Hermann and Oberguggenberger (see also ||l^) studied 
systems of singular, nonlinear ODEs in the Colombeau algebra. In this work we are going to use similar techniques 
to treat the geodesic and geodesic deviation equation for impulsive pp-waves in the Colombeau algebra. Despite the 
nonlinearities involved in these equations (which in principle could lead to trapping, blow-up or reflection of solutions 
at the shock, cf. |l^]) we are able to prove existence and uniqueness of geodesies crossing the shock hypersurface. We 
derive the (regularization independent) distributional limits of these solutions, making use of the notion of association 
(see Sec. H below) in the algebra, thereby significantly generalizing the results of In particular, the regularization 
of the (5-like wave profile will no longer be restricted to a "model ^-net" but belong to the largest "reasonable" class 
(cf. Definition |l] below). Moreover, note that the regularization idependence of the results has the following impor- 
tant physical consequence: in the impulsive limit the geodesies are totally independent of the particular shape of the 
sandwich wave. Hence the impulsive wave "totally forgets its seed" (cf. also the results in JT9|). 

Finally, we discuss the case of a nonsmooth wave profile / and give an outlook to current research which allows to fit 
our previous calculations into a manifestly covariant concept of Colomebeau algebras on manifolds. 



II. MATHEMATICAL FRAMEWORK 



A framework that allows consistent treatment of nonlinear operations with distributions and at the same time offers 
a well-developed theory of (linear and nonlinear) partial differential equations is provided by Colombeau's theory of 
algebras of generalized functions (cf. e.g. |p^- p^ , po[ ) . To begin with, we give a short description of the algebra we are 
going to use in the sequel. Let 



A(R") ^{ipe P(K") : J ip{x)dx = 1} 

Aq{W') - {(^ e ^o(K") : / 'p{x}x°' dx^O,l<\a\<q} {q E 



and set (for any fl C R" open) 

£{n) = {R : A(K") xn^C : x^ Ri<P,x) £ C°°(f})V(^ G -4o(K")} 

£m{^) = {ue £{n) -.^K ccnVaE 3p G No e Ap{m.") 

3c> 3r] > sup \d°'u{(p^,x)\ < ce^P (0 < e < ?/)} 
xeK 

j\f{n) ^{ue £{n) : V/C CC 17 Va e NJ^ 3p e No 37 e r Vg > p 

Vip e AgiR"") 3c> 3r; > sup \d°'u{ip„x)\ < ce''^'?)-^ {0 < e < t])} , 

where F = {7 : No ^ K+ : 7 strictly increasing, lim j{n) — 00}. Derivation d" is carried out with respect to x, 

n — *oo 

while the ip are treated as parameters. Also, for ip e T>{M."), (pe{x) = £~"(p(|). Note that (p^ ^ 6 in T>'{M."-). 
Elements of £m{^) are called of moderate growth. With pointwise operations £m{^) is a differential algebra and A/'(r2) 
is an ideal in £m{^)- The quotient algebra 

g{n) - £M{n)/M{n) 

is called the Colombeau algebra over C R". Elements of g{il) will be denoted by i? = c\[{R{lp, .))ip(zj,g] where 
{R{p, ■ ))(peXo is an arbitrary representative of R (again emphasizing the fact that the (^'s are viewed as parameters). 
For 17 = R" the map 

L : £'{fi) g{n) 

w d[{w * (p)^eAo] 
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(where * denotes convolution) is a linear embedding commuting with partial derivatives and coinciding with the 
identical embedding / — ^ cl[(/);pg^g] on 

Q is a fine sheaf of differential algebras on M" and there is a unique sheaf niorphisni T : V ^ Q coinciding with t 
on every £'{il) and rendering C°°{^1) a faithful subalgebra of G{^1). From the definitions it is clear that any element 
of is uniquely determined by the values of any representative on (/j^ for if g Ap with p arbitrarily large and e 
arbitrarily small (i.e. by its 'germ'), a fact that turns out to be very helpful e.g. in constructing solutions to differential 
equations in Q. 

Inserting points into elements of t/(M") gives elements of the ring of generalized numbers C(n), defined as C{n) = 
£{n)/M{n), where 

£{n) = {u : A(K") ^ C : 3p £ No G ylp(M") 3c> 3?? > 

\u{lp^)\ < ce-P (0 < £ < 77)} 
M{n) = {u : A(K") ^ C : 3p e No 37 e r Vg > p e 
3c> 377 > \u{ip,)\ < ce''('?)-f (0 < e < 77)} 

Thus elements of Q{W^) take values in C(77,). Explicit dependence of the ring of constants on n can be avoided by a 
more refined construction of the sets Aq in the definition of Q (see ) . Clearly, C ^ C via the canonical embedding 
c^cl[(c)^e^J. 

Componentwise insertion oi R ^ Q into a smooth function / yields a well defined element f{R) of if / is slowly 
increasing, i.e. if all derivatives of / are polynomially bounded. Moreover, if R is locally hounded, i.e. if it possesses 
a representative such that R{ipe, . ) is bounded uniformly in e on compact sets (for ip £ Ap{M"), p large) then f o R 
exists for any smooth /. 

Finally, we mention the notion of association in R2 G G{^) are called associated to each other (i?i w R2) 

if there exists some p G N such that Ri{(pe, ■ ) — -R2(v?e, . ) ^ in 'D'{D,) as e ^ for all tp G ^p(R"). In particular, 
if i?2 G I?'(f2) then R2 is called the macroscopic aspect (or distributional shadow) of Equality in V is reflected 
as equality in the sense of association in while equality in is a stricter concept (for example, all powers of the 
Heaviside function are distinct in the Colombeau algebra although they are associated with each other). 



III. EXACT SOLUTIONS OF GEODESIC AND GEODESIC DEVIATION EQUATIONS 

As in fl^ we consider the impulsive pp-wave metric 

ds^ = fix') d{u) du^ -dudv + {dx')^ , (2) 

where / is a smooth function of the transverse coordinates x' {i = 1,2). Our aim is to derive solutions to the 
corresponding geodesic and geodesic deviation equations in the Colombeau algebra. 

The general strategy for solving differential equations in Q is to embed singularities (in our case: 5) into Q which 
amounts to a regularization and then solve the corresponding regularized equations. In order to obtain general re- 
sults we are therefore interested in imposing as few restrictions as possible on the regularization of 5. The largest 
"reasonable" class of smooth^ regularizations of 5 is given by nets (pe)eg(o.i) of smooth functions satisfying: 

(a) supp(p,) ^ {0} (£ ^ 0) , 

(b) / Pe{x) — > 1 (e — > 0) and 



(c) 3t]>03C>0: J \pe{x)\dx < C Ve G (0,77) 



(cf. the definition of strict delta nets in ||l6|, ch. 2.7). Obviously any such net converges to S in distributions as e ^ 0. 
To simplify notations it is often convenient to replace (a) by 

(a') supp(/9e) C [-e,e] Ve G (0, 1). 

This motivates the following (cf. ITJ 



'^Note that since V is dense in practically even discontinuous regularizations (eg. boxes) are included. 
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Definition 1 A generalized delta function is an element D o/t/(R") possessing a representative {D{ip, . ))ipi£Ao such 
that 3p e No G Ap{W') 37] = r]{tp) > : 

(i) supp(i?(<^e, .)) C [-£,£] Vee (0,77) 

(ii) / D{ip^,x)dx 1 (e ^ 0) 



(iii) 3C = C{ip) > such that / |i:'(<^£, a;)| dx < C Ve £ (0, r/) 



The canonical embedding R = l{S) of course falls into this class but clearly there are many generalized delta functions 
that do not correspond to any distribution via l. Moreover, every generalized delta function is associated to S, i.e. all 
generalized delta functions equal S on the distributional level. In a sense, they may be viewed as 'delta distributions 
with a more refined microstructure' (fixing the additional nonlinear properties of the singularity). 
Again, condition (i) in definition |l] has been chosen in order to avoid technicalities in the proofs of the following 
results, which, however, remain true if (i) is replaced by 

(i') supp(i?(^e, .))^{0} (e^O) 
Finally, we need the following technical preparation (which is actually a generalization of appendix A of p^). 

Lemma 1 Let g : M" — > M", h : R M" smooth and let {pe)ee{o.i) be a net of smooth functions satisfying (a') and 
(c). For any xq, xq G R" and any e G (0, 1) consider the system 

Xe{t) = g{x,{t))pe{t) + h{t) 
Xe(-l)^Xo (3) 
Xe(-l) = Xq 

Let b > 0, M — J^^ J^^ \h{r) \ drds, / = {x G R" : |a; — XqI < + jiol + M} and a = mui ^ +\±^\^ 2LC 

with L a Lipschitz constant for g on L. Then ^) has a unique solution on = [— l,a — e]. Consequently, for s 
sufficiently small x^ is globally defined and both x^ and x^ are bounded, uniformly in e, on compact sets. 

Proof. The operator / Af , 

t s t s 

Af{t)^xo+xo{t+l) + J J g{f{r))p,{r)drds + J J hir) drds 

-1-1 -1-1 

is a contraction on the complete metric space 

{/gC(J„R") : \f{t)-xo\ <b + M+\io\} 

□ 

Let us now turn to the geodesic equation for the pp-wave metric (^. Using u as an affine parameter (which excludes 
trivial geodesies parallel to the shock) we obtain (cf. [p^): 

v{u) = f{x'{u))6{u) + 2d, f{x\u))x' {u)5{u) 
i\u)^]^d,f{x\u))6{u) 

Since all operations appearing in (^ are well-defined in Q (cf. the remarks following Theorem below), we may seek 
solutions of the corresponding initial value problem in the Colombeau algebra by embedding 5{u) into Q. In fact, it 
turns out that for any generalized delta function there exists a unique solution. Denoting the generalized functions 
corresponding to x^ and v by capital letters we state the following 

Theorem 1 Let D G Q(M.) be a generalized delta function, f G C°°(R^) and let vq, vq, Xq, Xq G R (i = 1,2). The 
initial value problem 
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V{u) =f{X\u))b{u) + 2dJ{X\u))X'{u)D{u) 

r{u)^^dJ(X\u))D{u) 

V{-l) = vo XH-1) = 
l/(-l)=«o ^'(-1) = il 

has a unique locally hounded solution {V, X^ , X"^) G C/(R)^. 

Note that we impose initial conditions in u = — 1, i.e. "long before" the shock. Choosing initial conditions at u = 
would mean to start "at the shock," which inevitably leads to regularization dependent weak limits. 
Proof. Existence: Choose p G N as in definition 0, fix G and let e < ?7(</j). Then componentwise we obtain 

the equations 

V{^,,u) ^ f{X\^,,u))b{^,,u) + 2d,fiX\cp„u))X' (cp^u) D{cp,,u) 

X'{pe,u) = ^dJ{X\p„ u)) D{p„u) 

V{ip,,-l)=Vo X\^,,-l) = xl, 
V{p,,-l)=vo X\^„-l)^xl 

According to Lemma |l|, the second line of (^) has a unique globally defined solution X^{ip^, . ) with the specified ini- 
tial values. Inserting this into the first line and integrating we also obtain a solution V{ip^, . ). From the boundedness 
properties of X'^{(pg, . ) established in Lemma |l| and the fact that {D{(p, . ))if,£Ao G Sm{^) it follows easily by induction 
that {X'^{(p^, . ))^£Ao ^-nd {V{ip^, . ))ip£Ao S'l'G moderate as well. Hence their respective classes in Q{M.) define solutions 
to (I). 

Uniqueness: Suppose that Vi = c\[{Vi{(pe, .))ip^Ao] ^-nd XI = c\[{X^{ipe, .))^i£Ao] are locally bounded solu- 
tions of (||) as well. On the level of representatives this means that there exist M — c\[{M{(pg, .))^p(zAo], A^' = 
cl[(7V'((^j, . ))^!zAo] G A/'(M) and n^^, n^,,, Uy, Uy G Mil) with 

i>i(^„u) = f{X{{^„u))D{^„u) + 2d,f{X{{^,,u))Xl{^,,u) 
D{ip,,u) + M{ip,,u) 

Xii^s,u) ^^dJ{Xl{ip„u))D{^,,u) + N\ip„u) (7) 
Vii'Pe, -1) ^vo+ ny{ipe) Xl{ip^, -1) = + n^^{ip^) 

Vl((^£, -1) = Wo + nv{^s) X{{'fe, -1) = ±0 + "i'lV'e) 

We have to show that {{V -Vi){(pe, ■ ))^eAo and {{X^ ~Xl){Lp^, . )),^eAa belong to the ideal J\f[M). Since G A/'(R) 
it follows that for p sufficiently large, e small and (f G ^p(M), N'^{(p^, . ) is bounded on compact sets, uniformly in e. 
Thus by Lemma ^ the same holds true for XI {ip^ , . ) and its first derivative. From we conclude 



US US 

+ \j j D{^„r)[d,f{X\^,,r)) - dd{Xl{^„r))]drds - J J N{^,,r)drds 
-1-1 -1-1 

Hence VT > 3p G Nq ^7 G F Vq > p V(y3 G y^,(M) 3C > 3?? > Ve G (0, 77) Vu G [-T, T]: 

u u 1 

\{X^-Xl){^,,u)\<Ce''^^^-^ + \j j JmfiaX\p,,r)) + 

-1 -r 

+ {l-a)Xl{ip„r))\da \ (A' - XI ) {ip, , r) | \D{ip„ r) \ dsdr 

By the boundedness properties of A* and XI and by (iii) , an application of Gronwall's Lemma to the above inequality 
yields the A/'-estimates of order for (A' ~ XI). A similar argument applies to the first derivatives. The estimates of 
higher order then follow inductively from the differential equation, so ((A* — Xl){(p^, . ))(p£Ao ^ A/'(R). Inserting this 
into the integral equation for {V — Vi), the A/'-estimates for {V — Vi) also follow inductively. □ 
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In the proof of Theorem |l| we have only made use of properties (i) and (iii) of the generalized delta fmiction D. On 
the other hand, property (ii) will be essential for the explicit calculation of distributional limits of the unique solution 
constructed in Theorem n^, cf. Sec. IV. Also, note that we did not have to impose any growth restrictions on / to 



obtain a well-defined element f{X^) of Q. This is of course due to the fact that any componentwise solution of the 
initial value problem necessarily is bounded, uniformly in e, on compact sets (for e small). 

Our next goal is an analysis of the Jacobi equation for impulsive pp-waves in the framework of algebras of generalized 
functions. As in to keep formulas more transparent we make some simplifying assumptions concerning geometry 
(namely axisymmetry) and initial conditions. Writing x — and y = x^ we suppose that / depends exclusively on 
the two-radius -y/a^^ + y^ and work within the hypersurface y = (corresponding to initial conditions yo = = yo). 
Furthermore we demand = = .tq. As was shown in in this situation the Jacobi equation 

= -Rl,T''T''N- , 

where N°-{u) — (N'^{u), N"" (u), N'^{u), (u)) denotes the deviation vector field takes the form 
iV" = 2[N''f'{x)S] ■ - Ny{x)S + [N''f{x)S]-~ N''f"{x)±^6 - 7V"/'(x)£(5 

N- = [N-f{x) + 1 N'=rix)]6 + i fix)N^S (9) 
Ny ^N"" = , 

where x is determined by (jj). Existence and uniqueness of solutions to the corresponding initial value problem in 
the Colombeau algebra is established in the following result where, for the sake of brevity we denote the ^/-functions 
corresponding to N'^ again by N'^. 

Theorem 2 Let D G Q{M.) be a generalized delta function, f S C°°(M), n° G and let X denote the (unique) 
solution to system ^) with initial conditions and simplifications as discussed above. The initial value problem 

TV" = 2[N''f'{X)D]--N='f'{X)D+[N''f{X)Dy- 
-N''f"{X)X^D - N'^f'{X)XD 

= [N-nX) + i N-f"{X)]D + i f{X)N^iD ^^^^ 

NV ^ TV" = 
7V"(-1) = n'^ iV°(-l) = ii" 



has a unique solution N"" G Q{M.)'^. 

Proof. Since the equations are linear in the components of the deviation field we are provided with globally defined 
solutions on the level of representatives. The last two equations are actually trivial and so is the first one once we 
know that its right hand side belongs to t/(M). Hence we are left with the equation for which is of the form 
N{t) = f"{X{t))D(t)N{t) + H{t) with H in ^(M). Using the boundedness properties of X established in Lemma |l| 
the ^M-bounds for easily follow from Gronwall's lemma. 

Uniqueness is established along the same lines again using Gronwall-type arguments. □ 
In the above proof we have again only used properties (i) and (iii) of the generalized delta function D. 
To conclude this section we remark that unique solvability of the geodesic and geodesic deviation equation for (^) is 
not confined to the case where the profile function / is smooth. Indeed, it turns out that for a large class of generalized 
profile functions (those that are not "too singular" ) Theorems Q and || retain their validity. More precisely, we have 
to demand that / belongs to the algebra of tempered generalized functions ]l5[ | to make sure that the composition 
f(X) is well defined and that VV/ is of L°°-log-type | p7|jl8| to ensure existence and uniqueness of solutions to ^ 
and (|l^). However, to include many physically interesting examples (cf. ^2l|]') one has to cut out the worldlinc of the 
ultrarelativistic particle, i.e. the w-axis from the domain of definition (cf. |p2[). 



IV. DISTRIBUTIONAL LIMITS 



In this section we are going to calculate the distributional limits (or, in the terminology of Colombeau theory: 
the associated distributions) of the unique solutions to the geodesic and geodesic deviation equation constructed in 
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Theorems |l| and ||. In |jT^ distributional limits for regularized versions of these equations have been calculated using 
a model delta net regularization. Translated into our current setting this amounts to using the particular generalized 
delta function D = t((5). Our aim is to extend the validity of the limit relations derived there to the case of solutions 
in the Colombeau algebra and to generalized delta functions. At the same time we will be able to prove stronger 
convergence results in some cases. 

Theorem 3 The unique solution {V,X'^) oj the geodesic equation satisfies the following association relations: 

X'^xl^xl{l + u) + ]^d,f{xl + xl)u+ (11) 
V ^v^ + ^o(l +u) + f{xl + xl)9{u) + d,f{xl + ±1) (^4 + \d,f{xl + 4) j u+ (12) 

In addition, if = c\[{X^{ip, . ))ipeAa] then 3p G Nq such that \/(p G Ap 

X\ip,, u)^xl,+ i^(l + u) + ^dJixl + xl,)u+ (13) 
for e — !■ 0, uniformly on compact subsets of K. 

Proof. Choose p g No as in definition |l| for D, let if £ Ap and e < r]{ip). Since integrating amounts to convolution 
with the Heaviside function, which is a continuous operation on the convolution algebra of distributions supported in 
a cone, in order to prove ( p^ ) it suffices to show that 

X\Ve, .) = ^d,f{X\^,, .))Di^,, .) ^ ldaixl + il)S 

in distributions. We first note that X{(p^,et) — > Xq + Xq uniformly as can be seen from the integral equation for X* 
(cf. (11) in Now iiij e V{R) then 

m^^f{X\^,, t))D{^,,t)dt ~ dj{xl + i^)^(O) I 



< sup I V(i)a^/(XX¥'e,t)) -5./(x^ +i^)V(0) I / I D{ip,,t) I dt + 

-e<t<e J 

— e 

e 

+ J I D{ip,,t)dt - 1 I dJix'o + ij,)^(0) 

— e 

So the claim follows from properties (iii) and (ii) of the generalized delta function D. Since X'^{(pe, t) is bounded on 
compact sets, uniformly in e, it follows that the family {X'((/3£, t) : e £ (0, 1)} is locally equicontinuous. Hence Ascoli's 
Theorem implies (|l3[). Concerning (p^), as above it suffices to calculate the limit of 

V(^,,u) = [f{X\^,,u))D(^,,u)] + d,f{X\^,,u))X' {ip„u)D{ip,,u) 

whose first summand converges to /{x'q + Xq)S by an argument similar to the one above. For the second summand we 
have 

dJiX\^„u))X\ip„u)D{^,,u) = dJiX\ip,,u))D{^,,u)xl, 

" v ' 

(*) 

t 

+ ^d,f{X\ip,,u))D{^,,u) J dJ{X\^„s))D{^„s)ds 

— e 

and (*) dif{xl + xI)xq5. Finally, since 
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-i5./K + ij,)V(0) / D{^,,t)dt-.0, 

the claim follows. □ 

In calculating distributional limits for 

the solution of the Jacobi equation to maintain clarity of formulae we shall make simplifying assumptions on the 

Ar°(-1) = (0,0,0,0) , ^ 

initial conditions, i.e. . (14) 

iV"(-l) = (a,6,0,0) . 

Then we have 

Theorem 4 The unique solution of the geodesic deviation equation ^lul) satisfies the following association relations: 



iV-«ia/'(xo)K + 0(u)) (15) 
« b{l + u) + a[f{xo)6{u) + ^f'{xof{d{u) + u+)] (16) 

Proof. The general structure of this proof is 'isomorphic' to the calculation of distributional limits for the regularized 



Jacobi equation in |12|. The main difference is that for representatives of generalized delta functions dominated 
convergence arguments are not applicable which makes the calculations more tedious. Nevertheless, using the uniform 
convergence of X^{(p^, .) established above, all steps carried out in |l^ can be adapted to the present situation as 
demonstrated in the proof of Theorem ^. □ 



V. DISCUSSION AND OUTLOOK 

In the previous section we have shown that the unique solutions to the geodesic and geodesic deviation equation 
in the Colombeau algebra possess a physically reasonable macroscopic (i.e. distributional) aspect: even within the 
natural maximal class of delta-regularizations (namely the class of all generalized delta functions) the regularity of 
the equations is sufficiently high to ensure distributional limits corresponding to physical expectations. More precisely, 
from the distributional point of view, the geodesies correspond to refracted, broken straight lines as suggested by the 
form of the metric. The scale of the jump and kink is given by the values of / and its first derivatives at the shock 
hypersurface, which can be traced back to the values at the initial point (u = —1), thereby precisely reproducing 
Penrose's junction conditions [Q. The distributional limit of the Jacobi field suffers a kink and jump in x-direction 
as well as an additional ^-pulse in w-direction, which may be understood from the form of the geodesies. For a more 
detailed discussion see ||l^ . 

Finally we make some comments on diffeomorphism invariance of our results. Whereas the fine sheaf of Colombeau 
algebras can be lifted to manifolds in a straightforward manner, the action of a diffeomorphism does not commute 
with the canonical embedding V ^ Q. The reason for this is that convolution relies on the additive group structure 
of R" and is therefore not invariant under the action of diffeomorphisms. Note, however that our calculations did not 
use the embedding via convolution and therefore are not affected by this defect. 

A solution to the above mentioned problem was proposed in using a modified definition of the moUifier spaces Aq. 
A key ingredient of this construction is that diffeomorphisms act on the ip^s, introducing an implicit a;-dependence 
into the first slot of the Colombeau functions R{(p,x). Hence, to retain smooth x-dependence of R the concept of 
Silva-differentiability was used. Future work will be concerned with a simplified concept of Colombeau algebras on 
manifolds using calculus in convenient vector spaces [ p4| . A main goal of this line of research is to provide a workable 
solution concept for singular differential equations on manifolds. 
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